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SUMMARY 


Three advanced methods to compute the rotor thickness noise which is pre- 
dominant in the case of high speed rotor have been developed. These methods are 
deduced from a previous method by transforming the integral coordinate, commut- 
ing the order of integration and differentiation, and/or performing chordwise 
integration analytically with some adequate assumptions. The necessary compu- 
tational times and waveforms obtained by the previous and three advanced 
methods were compared. It was then concluded that the advanced methods could 
save the computational time very much compared with the previous method in 
keeping the same accuracy. 


INTRODUCTION 


Farassat has proposed a method to calculate the thickness noise which is 
produced by moving bodies of finite volume or thickness normal to the moving 
direction and applied it to the rotor noise of helicopters (ref.l). 

Based on his work, the authors made clear the cause and the characteris- 
tics of the rotor rotational noise. Through these studies (ref. 2 and 3) with 
the concept of the "influential surface", which is an integral region at 
retarded time, and of the distribution of source strength, many acoustic char- 
acteristics of helicopter rotor noise have been clarified by the analytical pre- 
dictions and numerical calculations, and verified by the experimental tests. 

It has been shown that among the rotor rotational noise components the rotor 
thickness effect executes a dominant roll rather than the loading effects such 
as thrust noise and drag noise in the case of high blade-tip speed. 

In the computational process of the thickness noise, however, the calcula- 
tion has the worst converging characteristics and thus needs a lot of computer 
time in getting solutions within enough accuracy. It has, therefore, been 
expected to develop an improved method of calculation for obtaining an advanced 
form of solution which is more convenient to get the result with less computer 
time and enough accuracy. 
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SYMBOLS 


B 

CO 

C h 

f=0 


g=o 


h 

K 

M 

r 

P 

R 

Ro 

r 

ro 


t 

V 

V 

n 

V 

n 

x 

y 


number of "blades 
sound speed, m/see 
chord length 

equation of body surface 

j r) 3 -h(m, n 2 ) at "the upper surface 

I -ri 3 -h(rii, 02 ) at the lower surface 

equation of acoustic sphere = T-t+[r] /cq 
blade thickness, m 

integrand of a modified solution of thickness noise, see equation 

(10) 


relative Mach number 
acoustic pressure, kg/m 2 
rotor radius, m 
blade cut-off radius, m 

distance between source and observer, =|x-y|, m 
distance between hub center and observer, m 

unit vector in direction of propagation in blade fixed coordinate 
f =(r , r , r ) T = (3r/3r) 1 , 3r/3q 2 , 3r/3n 3 ) T 

I n 1 b 2 b 3 

\ r 1 sin^-r 2 cosi/i r 3 

U< 7 — >°*-> 

observer time 

T 

rotor hub velocity, =(Vi ,V 2 ,V 3 ) , m/sec 
= ( V ,V ,V ) T , m/sec 

n 1 n 2 B3 

blade element velocity normal to blade surface, m/sec 

T 

observer position vector, =(xi,x 2 ,X 3 ) } m 

T 

source position vector, =(yi,y 2 »y3) , m 
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I 


a 


n 

e 


p 

T 


angle between radiating direction and rotor plane, rad or deg 

curve of intersection of body with acoustic sphere 

position vector in blade fixed Cartesian coordinate system 

angle between radiating direction and normal direction to blade 
surface, rad or deg 

density of undisturbed fluid, kg*see/m 4 
source time, sec 


ip blade azimuth angle, rad or deg 

ft rotor rotational speed, rad/sec 

A directional parameter in modified method, see equation (8) 

X directional parameter in previous method, see equation ( 3 ) 


Superscript 



transposed of 


) 


Subscript 

[ ] . value at a retarded time 

ret 

Operator 

V gradient operator in fluid fixed coordinate, =9/9y^, l/m 


PREVIOUS ANALYSIS — METHOD (A) 


According to Farassat (ref.l), the rotor thickness noise, p, for a given 
observer position, X, and time, t, is given by 

p c 0 3 ( 3h/ 3n 1 )Vn 1 

p (x,t) = 2 V atl-Tilr r\ drdT ^ 


where h(ni) is the blade thickness distribution as a function of chordwise 
coordinate, n 1 - Here the relative speed of a blade element with respect to the 
fluid and the directional parameter. A, are respectively given by 


V 

m 


Visini|j-V2Cos^+ri2^ 5 


( 2 ) 
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A = |Vf|sin0 = |(-3h/3n 1 , -3h/30 2 , ±l)^jsin0 

( 3 ) 

= {1-r 2 +( 3h/ 3 n i ) 2 (l-r 2 
03 01 

where r is the distance between the observer and the source, and Vf means the 
vector outward normal to the blade surface. Other symbols related to the above 
equations are listed in SYMBOLS. 

The integrations appearing in equation (l) may be understood as follows: 

In order to calculate the noise field of a moving source, the concept of re- 
tarded time designated by [ ] re ^_ must be considered, because the pressure 
change generated at different points and times might be received by an observer 
at a given time simultaneously. It is, thus, useful to know all regions of 
noise sources that have influence on the observer at the given time. These 
regions will be formed by the loci or trajectories of closed curves called 
"T-curves" for the given observer time and position with the change of the 
source time, -°°<x<t . 


An external surface of each region has been named "influential surface" . 
All sources distributed on these surfaces must be integrated to give the in- 
stantaneous pressure of the given observer time and position. By considering 
the shape of the respective influential surface, many typical effects of 
various rotor operating parameters on each noise component have come to be 
predictable analytically (ref .2 and 3) . 

T 2 

The integration /^dlcgdi: should be performed on the influential sur- 
faces which are the loci of intersections between an "acoustic sphere," 
g=t-T-[r ] re ^/cg=0, and blade surfaces, f(y,x)=0, for a given observer time, t. 

The times specified by Tj and T 2 are source times at which the acoustic sphere 
enters and leaves the blade respectively. With the lapse of the source time, 
the acoustic sphere contracts toward its center, just where the observer 
locates, with the speed of sound, eg , while the blade rotates around the rotor 
axis. Fig.l shows the geometric arrangements of an influential point. The in- 
fluential point is given as the intersection between an acoustic line and a 
specified point, ( 01 , 112)5 on the i-th blade and is determined by the following 
equation: 

F' = cosa* (0 2 cos^-R)^/co+^ 

-(fit+277(i-l)/B-n 1 /n 2 ) = 0, 

where 


O<0i5C h 

R q ^ 02< r - 


(5) 


The angular velocity of the influential point, 3^/3t, which plays an im- 
portant role in the noise calculation as will be stated later, can be derived 
by differentiating equation ( h ) as 


3i|//3t = fi/(l-M^), 


( 6 ) 


326 



where M r is relative Mach number of the specified blade position with respect 
to an observer. Since M r has positive peak at ^=90°, 3tjj/3t takes maximum value 
at this azimuthal position. 

In the numerical calculation by the previous method, factors which decide 
the precision level of calculation are: i) mesh dimension of the Simpson sum, 
ii) dimension of the observer time increment, iii) precision rank adopted in 
the computation. After some trial computations by using a computer, FACOM 
230/75 > it was concluded that the double precision was necessary and sufficient 
for the present calculation and that the number of integrating point for ob- 
taining reliable results was about 10 4 , which corresponded to the following 
dimension for one element: The tangential partition was (chord length) /20; the 

radial partition was (rotor radius)/300; and An=0.l6. This mesh dimension re- 
quired about twenty seconds to compute the total pressure value of a given 
observer time and position. Then it was decided that the number of discrete 
observer times in the one blade passing period was thirty six in order to give 
a total computational time within twenty minutes. 


MODIFIED ANALYSES — METHODS (b) AND (C) 


Undesirable defects on the accuracy in the previous computation were caused 
by the numerical differentiation, 3/3t, and the numerical double integration on 
the influential surfaces, //dr dr. If the observer time differentiation, 3/3t, 
can be put inside the integration, and be performed analytically, then the 
numerical differentiation will disappear and very much computational time can 
be saved. 

In the previous analysis, the integration on the influential surface was 
performed along a contracting acoustic sphere, dldx. Here, in the present 
analyses, by using a polar coordinate system associated with the rotor disk, an 
elemental area of the integration, ri 2 d^dn 2 5 can be transformed as follows: 

codrdT n2<tydri2 

X/|vf J = A ^ 


where 


A = (l+M2-2M n cos9) 1/2 

M n = v n /c 0 (8) 

V n = -|f/l Vf l 

Then, by lasing the above relation, the solution of wave equation of the rotor 
thickness noise can be rewritten in the new coordinate system as follows: 
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Tl2 dr l2 


(9) 


where 


p(x * t> = -fe !QCr ] 

u ^ ret 


K = 3^ 1 Vni/r l Vf l A - (10) 

In the ranges of integrations, ipi and 1^2 are the azimuth angles of the leading 
edge and trailing edge of the influential surface respectively, which are 
functions of the span position, p , the observer time, t, and position, x; and 
R and Ro are the rotor radius and the blade cut-off radius respectively. The 
integrand, K, is considered to be a function of independent variable, m or i(j , 
n 2 , t, and x. Thus, the double integration, //n 2 ^dn 2 > gives an influential- 
surface integration of a single blade. 

Fig. 2 shows two different arrangements of integrating points on the same 
influential surface. It can be seen that the present methods have a well- 
fitting coordinate converging to the integral region. The method (B) is derived 
from equation ( 9 ) by performing numerically the differentiation and the 
integration. 

By commuting the order of the differentiation and the integration, equation 
( 9 ) becomes 


p(x,t) 


P 

2tt 



'h (t) 
* 2 (t) 


3K( ^ ) 

3t 


!d^ -K(^ 2 ) 
'^,T12 


d\p2 


d] P 1-j 

+K ^ 1 ^9tJ Tl2dn2 ’ 

ret 


where 9K(^)/9t should be considered as follows: 


(11) 


9K 


’ 3K 

+ M anj' 

at 

L,n 2 

. 9t 

9n 1 3t 

ni 


9K _9r 
9r 3t 


3K 

9A 


9A 

3t 


ni ,n2»>P 


7 9K 9 2 h 

_^9( 9h/9ni )9nT 


( 12 ) 


9K 9Vn ! | 3K 9r 9K 9 j Vf ■[ ( 3K 9A ] 9m 

9V 9ni 3m 9 J V f [ 9ni 9A 3nV 3t 

111 n 2 ,<P * 


The method (C) is derived from equation (ll) by performing the numerical 
integration . 
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ANALYTIC METHOD METHOD (D) 


Consider the following case for simplicity of analysis: 


(i) hovering state. 



= n 2 fi 


and 


(13) 


(ii) far field, R<<r that is r=rg=const., A=l, and Vf|=l. 


Then the kernel given by equation (10) and the acoustic pressure can be 
written respectively as follows: 


K = ( 3h/3n i )n 2 ^/ r 0 


(lb) 


and 


p(x,t) 


2irr 0 



dn2 

ret 


2itr 0 

2Trr 0 


R 

Ro 

R 


n: 


Vl 


V2 

<h 


( ( 

Ro ^2 


( 

L 3t 3n 1 

3 2 h 


dV+ 


3h 


3n i 


J V 

3 2 h 3n l 


3Vl 

It 


VI 


3t 3n i 3n^ 3t 


r 


Mv+ 


3h 

3ni 

3h 


3ni 


3^2"! 

It 

V2 

3Vi 

Vi3t 


dn2 

ret 

3h 


(15) 


3n ■ 


3^21 
V 2 9t 


dri2 ■ 

ret 


In a double parabolic airfoil, for example, the thickness change can be 
given by 


3h _ 
3n 1 


2h 


max „ ,„ \ 

~ d-2m/c h ). 


which is, at the leading edge, n^O OJH^)* and the trailing edge, 
n i =Ch(V=V 2 ) , given by 


3h 
3n 1 



3h 
3n 1 



2hjnax 



(16) 


(IT) 
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Combining equations (l6) and (4) yields 


_ 2h ma x r •, 

am' c, L x ~ 


2H2 


2tt, 


{ftt+ f-(i-l) 

C h B 


(18) 


-ip- — ( t ^ p 03 ^”^) 2030 ^ 3 

c o 


By substituting the above two relations into equation (15), the acoustic 
pressure can be obtained in a decomposed form as follows: 


p(x,t) = pi(x,t)+p 2 (X,t)+p3(X,t) 


where 


Pl(x,t) 

p 2 (x,t) 

p 3 (x,t) 


P^hmax 
’ r ° C h 'E 0 
t , R 

pfth 
ffr 0 C. 


1 *- 

; Ro 

( 

/ Tin 


2r\ 2^ 

— — 1-^2 ) >a.ri 2 

°h 

dH2 


(19) 


h ' R 0 
P_^h max I n 

' "D - 


irroC 


ri^dn; 


h R 0 


As written above, the chordwise integration and the time differentiation 
have been performed analitically . Three components of equation (19) are cor- 
responding to the first, second, and third terms of equation (ll) exactly. An 
example of computed waveform is shown in Fig. 3, in which solid lines show 
components of the respective term and a dotted line shows the total acoustic 
waveform obtained by summing up these three components. 

As shown in Fig.i+, the influential surfaces of multi-bladed rotor can be 
made from those of one blade at appropriate observer time. Then, once the sound 
pressure of one blade is computed, the acoustic pressure of the multi-bladed 
rotor can be easily obtained by summing up those of each blade one after 
another. Fig. 5 shows an example of this process for three bladed rotor. 


COMPUTATIONAL RESULTS 


The dimensions of an exemplified rotor are given in Table I. The computa- 
tional procedure is shown in Fig. 6. Shown in Fig. 7 is an example of the 
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computed waveforms "by four different methods, (A) through (D) . It can "be seen 
that these four methods give the very good coincidence in the waveform and the 
negative peak value, hut a slight difference in the positive peak value. The 
difference between the waveforms of (B) and (C) is too small to be distinguished. 
It is interesting to note that the area made by the positive acoustic pressure 
and zero line in one period is equal to that made by the negative one. 

Fig. 8 shows the change of computed pressure amplitude in methods (B) and 

(C) in relation to the spanwise mesh size for a given chordwise mesh size by 
circles, and to the chordwise mesh size for a given spanwise mesh size by tri- 
angles. The difference of sensitivity of the accuracy for the number of span- 
wise and chordwise partitions are obvious. The insensitive tendency of the 
chordwise partition may result mostly from the blade contour of symmetric profile 
without singularity. 

Fig. 9 shows the computed pressure amplitudes or peak values of pressure 
versus the number of spanwise partitions in four methods. As the number of 
partitions increases the pressure converges to an expected true value whereas 
the computation time increases predominantly as shown in Fig. 10. It can be 
seen that choice of proper coordinate, analytic chordwise integration, and the 
consequently simplified program used in the present methods saved the computa- 
tional time significantly. This guarantees that the present methods (c) and 

(D) need respectively only one tenth and one hundredth of the computational 

time of the method (A) to get the acoustic pressure in the same accuracy. 


CONCLUSION 


Four different methods to compute the rotor thickness noise have been 
compared. They are (A) the previous method based on equation (l); (B) the 

one modified method based on equation (9) in which the surface integration is 
performed in the rotor- fixed-coordinate system; (c) the other modified method 
based on equation (11) in which the numerical differentiation is discarded; 

(D) the still other modified method based on equation (19) in which the numeri- 
cal chordwise integration is further discarded. All methods have shown the 
good coincidence in both the waveform and the peak amplitude. 

The introduction of the rotor-fixed-polar-coordinate system has brought 
better characteristics in convergence and accuracy of the numerical integration 
than the previous method (A) in which the observer-fixed-polar-coordinate sys- 
tem was adopted. Specifically, the method (d), in which the chordwise integra- 
tion was performed analytically by assuming the hovering state of rotor and the 
far field location of observer, has saved further computational time. The 
computational time in the method (d) was only one hundredth of method (A) and 
one tenth of methods (b) and (c) in keeping the same accuracy. It can, thus, 
be concluded that the computing method of the rotor thickness noise was 
improved very much by introducing the advanced methods (.C) and (d) in the 
accuracy and speed. 
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TABLE I 


Dimensions : 


Rotor radius; R, m • • • 5.0 

Blade cut off; Ro ,m 0.7 

No. of blades; B 2 

Blade chord ; Ch, m 0.4 

Blade thickness ratio; h max /Ch • 0.10 

Airfoil section NACA 0010 

Operating conditions : 

Forward speed; V, m/sec 0 (hovering) 

Tip Mach number; Mt 0.9 

Observer : 

Observer position; ro ,m 50 

Observer elevation angle ;a^, deg 0 
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Figure 5.- An example to obtain waveform of multi-bladed rotor from 
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Figure 6.- Flowchart of modified mehtods to compute the rotor thickness noise. 
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Figure 7.- Waveforms given by four methods. 



0 5 10 20 40 80 160 

NUMBER OF PARTITIONS 


Figure 8.- Converging tendency with the number of partitions 
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Figure 9.- Converging tendency of the four methods. 
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Figure 10.- A comparison of computational time between four methods. 




